Scattering of atomic dark-bright solitons from narrow impurities 
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In this work, we study the collision of an atomic dark-bright soliton, in a two-component Bose- 
Einstein condensate, with a Gaussian barrier or well. First, we present the results of an experiment, 
illustrating the classical particle phenomenology (transmission or reflection) in the case of an equal 
barrier or well in both components. Then, motivated by the experimental observations, we perform 
systematic simulations considering not only the case of equal heights, but also the considerably more 
complex setting, where the potential affects only one of the two components. We systematically 
classify the ensuing cases within a two-parameter diagram of barrier amplitudes in the two compo- 
nents, and provide intuitive explanations for the resulting observations, as well as of their variations 
as the size of the barrier changes. 

I. INTRODUCTION 

Atomic Bose-Einstein condensates (BECs) [J [5] provide an ideal platform for the study of nonlinear phenomena 
at the mesoscopic scale (see, e.g., the reviews [3H5]). In this context, of particular interest are multi-component 
BECs, which - in their simplest version - are composed of two different hyperfine atomic levels of the same atom 
(e.g., ^^Rb [3 H]). In such systems, a rich variety of structures can be observed, which can not arise in single- 
component BECs. Regarding macroscopic nonlinear excitations of multi-component BECs, a distinctive feature of 
interest involves the potential formation of dark-bright (DB) solitons. This type of vector soliton consists of a dark 
soliton in one component coupled to a bright soliton in the second component. These solitons in repulsive BECs 
are usually referred to as "symbiotic" solitons, with this characterization stemming from the fact that the bright 
component cannot be supported in a stand-alone fashion (it is only supported as such in attractive BECs [S]), and is 
only sustained because of the presence of its dark-counterpart, which acts as an external trapping potential. 

Dark-bright solitons have been studied extensively in different settings in a large number of theoretical works (see, 
e.g., Refs. |5HT5]). while they have also been observed in experiments, both in two-component ^^Rb BECs [TM^ 
and in nonlinear optics [2^28) . In the BEC context, these experimental studies have chiefly involved the dynamics 
of a single DB soliton in a trap [HI [21] , the generation of multiple DB solitons in a counterflow experiment [20] , the 
study of their interactions [22], as well as the creation of SU(2)-rotated DB solitons, in the form of beating dark-dark 
solitons [13 m]. 

On the other hand, the fundamental problem of the interaction of solitons with localized impurities has been 
considered both in nonlinear wave theory [29^ and solid state physics [30 . The interaction of cither bright or dark 
solitons with (5-like impurities has been investigated in the framework of the nonlinear Schrodinger (NLS) equation 
(see, e.g., Refs. [5Tti55] ). Relevant studies in the physics of atomic BECs have also appeared some time ago (see, e.g., 
Refs. [5^H55] ). but also more recently, both in the setting of potential wells [3nilinij and in that of barriers [¥T1 142) : 
see also the very recent work of [43]. In this context, localized impurities can be created as focused far-detuned laser 
beams and have already been used in experiments with dark solitons [HJ US] ; we also note very recent experiments 
with matter-wave bright solitons of ^Li [15] and *^Rb [JT] atoms and localized barriers. However, such soliton-defect 
interactions are far less well explored in the case of the multi-component setting (see, e.g., Ref. [T7] where the statics 
of DB solitons was studied in the presence of 5-like impurities) . It is the aim of this work to address this problem and 
study, in particular, the scattering of atomic DB solitons at narrow impurities. 

Our presentation will be structured as follows. In section II, we present the relevant prototypical model setup in 
the form of two coupled Gross-Pitaevskii (CP) equations describing the dynamics of a binary BEC with repulsive 
interactions; DB solitons for this model are presented as well. We also present results of a prototypical experiment 
where the scattering of DB solitons at a barrier - which is present (and equal) in both atomic components - is 
studied; this experiment provides the motivation for a more systematic theoretical study which is presented in the 
next section. In particular, in sections III, we numerically explore the dynamics of single DB solitons in a harmonic 
trap. The simpler scenario that is studied refers to the case where the impurity is equal in the two components 
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- as in the case of the experiment; we find that a particle-based phenomenology is sufficient to capture the main 
characteristics here. On the other hand, we identify a far more significant wealth of possibilities in the setting where 
the barrier (or well) is applicable only in one of the two components. We present a systematic study within the 
plane of the amplitudes of the two-components, providing intuitive explanations (on the basis of effective potentials), 
wherever possible, for the observed phenomenology. Finally, in section IV, we summarize our findings and present 
our conclusions, as well as a number of directions for potential future studies. 



II. MODEL AND EXPERIMENTAL MOTIVATION 
A. Gross-Pitaevskii equations and dark-bright solitons 



We consider a two-component BEG composed of two different hyperfine states of the same alkali isotope. If this 
binary condensate is confined in a highly anisotropic trap (with longitudinal and transverse trapping frequencies 
ujx <C wj_), then the mean-field dynamics of the BEG can be described by the following system of two coupled CP 
equations pQ[2]: 

ihdt^Pj - + ■'^j(^) - + J29jM^^ (1) 

where (a;, t) (j = 1, 2) are the macroscopic wave functions of the two components normalized to the numbers of atoms 
Nj = l^jl'^dx, m is the atomic mass, are the chemical potentials, Qjk = 2hujj_aji~ are the effective ID coupling 
constants {ujk are the s-wave scattering lengths), while Vj{x) denote the external trapping potentials for each species. 
In our considerations below, we will assume that the component 1 (2) supports a dark (bright) soliton; additionally, 
we will assume that both components are confined by the usual harmonic trap, namely Vh{x) — (l/2)ma;^a;^, while - 
for each component - an additional localized "impurity" potential, which can be generated by off-resonant Gaussian 
laser beams, is also present. Thus, the external potentials Vj{x) for each of the two components are described as: 

^ 2x'^\ ( 'lx^\ 

yi{x) ^Vh + Eaexp ^-^j - ^2(2;) ^ Vh + ^^exp ^-^j , (2) 

where the parameters Ed, Eb and e^, £b set, respectively, the amplitudes and widths of the impurities in each compo- 
nent. Notice that for a blue- or red-detuned laser beam, the impurity potentials can either repel {E^fi > 0) or attract 
{Ed,b < 0) the atoms of the respective component of the condensate. 

We now cast Eqs. ([l]) into a dimensionless form as follows: measuring the densities length, time and energy 

in units of 2aii, a± — ^/hjuj±, ljJ^ and huj±, respectively, Eqs. ([ij are reduced to the form: 

idtU = -^dju + Vi{x)u + {\u\^ + gi2\v\^ - fii)u, (3) 

idtV = -^dlv + V2(x)v + {gi2\u\^ +922\v\^ - f^2)v. (4) 

In Eqs. ([3|-(|4|, the wave functions u and v correspond to and -02 respectively, the normalized nonlinearity coeffi- 
cients are given by gj2 = 5^2/311, while the normalized harmonic trap potential (incorporated in Vi{x) and V2{x) as 
discussed above) is now given by Vh{x) — (l/2)ri^a;^, where fl — uj^juj^. 

Notice that in the GP Eqs. ([3|-(|4| the number of atoms Nu^b of each component is conserved; in fact, Nd,b are 
given by N^^b = {o.^_/'2.aii)ND,Bi where Nd — |u(a;)pda; and Nb = |w(a:)pdx are the respective integrals 
of motion of the normalized GP Eqs. (|3])-(|4|. 

As mentioned above, in the physically relevant setting of '^^Rb, the scattering lengths characterizing the intra- and 
inter-component atomic collisions are almost equal; thus, to a first approximation, one may assume that gi2 = 322 ~ Ij 
which means that the system of Eqs. (|3|-([4| is of the Manakov type j48j; in this case, the system is integrable in 
the absence of the external potentials Vi,2(x) and admits exact analytical dark-bright soliton solutions. Particularly, 
considering the boundary conditions |up — >■ jj-i and jwp — )■ as \x\ — 00, Eqs. (|3|-(|4| possess an exact analytical 
one-DB-soliton solution of the following form (see, e.g., Ref. [S]): 



uj:i^{x,t) = y//ii{cos 0tanh^ + i sin (/)}, (5) 
VT>B{x,i) — ?7sech^exp[iA;a; -I- z6'(t)], (6) 
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FIG. 1: (a-d) Time evolution of the soliton-barrier interaction. The peak barrier potential is ^bec ■ The in-trap evolution 
times after generation of the solitons are a) 250 ms, b) 500 ms, c) 750 ms, d) 1000 ms. (e-h) Similar to (a-d) but for a peak 
barrier potential of 1.11 • fisEC and evolution times e) 300 ms, f) 600 ms, g) 700 ms, h) 1000 ms. In all cases the chemical 
potential of the BEG, hbec, is approximately 36 nK (about 750 Hz). 



where ^ — D[x — a;o(t)), (j) is the dark soliton's phase angle, cos0 and rj represent the amplitudes of the dark and 
bright solitons, and D and XQ{t) are associated with the inverse width and the center position of the DB soliton. 
Furthermore, k = Dtan0 = const and 6(t) are the wavenumber and phase of the bright soliton, respectively. The 



above parameters of the DB-soliton are connected through the equations: D — /ii cos 
9{t) = (1/2) (D^ — fc^ + /i2 — where xq is the DB soliton velocity. 



rj , Xf) = Dtancj), and 



B. Experimental results 



Having introduced our setup, we now proceed by presenting results of an experiment dealing with scattering of 
atomic DB solitons at barriers. In fact, the results that will be presented below, motivate the more detailed theoretical 
investigation of this paper, but also illustrate the experimental tractability of this direction, and showcase prototypical 
results along this vein. 

Our experimental results are summarized in Fig. [T] There, it is shown that, depending on the barrier height, DB 
solitons are either reflected by or transmitted through a repulsive barrier. The experiment is conducted with a BEG 
of 4.5 • 10^ atoms of *^Rb confined in an optical dipole trap with trapping frequencies {waxiai, i^vcrticai, (^horizontal} = 
2tt X {1.4, 120, 174} Hz. The solitons are generated by transferring a small fraction of the atoms from the initial 
\F,mp) = |1,— 1) hyperfine state to the |2,— 2) state and exploiting a counter-flow induced modulational instability 
[20| generated by a magnetic gradient along the axial direction. The number of solitons, as well as their initial 
positions, can be controlled by adjusting experimental parameters such as the number of atoms transferred into the 
second state and the strength and duration of the magnetic gradient used to induce the counter-flow. For the data 
presented in Fig. [T] the solitons are generated in the left part of the BEG, the magnetic gradient is subsequently 
turned off, and the solitons start moving towards the trap center. The oscillations of individual solitons in a trap have 
been investigated in detail in |21j . For the present data, we additionally ramp on a repulsive barrier at the center 
of the trap. The barrier is generated from a 660 nm laser beam with a narrow waist of approximately 18 /^m in the 
direction of the BEG axis and has an aspect ratio of 4. For imaging, the two components of the BEG are vertically 
separated during a short free expansion time of 7 ms for the upper cloud and 8 ms for the lower cloud [20', '21] . For 
barrier depths larger than the chemical potential [cf. Fig. [T|e-h)] we observe confinement of the dark-bright solitons 
to the left half of the BEG. This is consistent with having two isolated BECs. For a barrier depth of approximately 
half the chemical potential [cf. Fig. [ij a-d)] we observe solitons penetrating through the barrier; see, e.g., especially 
the panel (c) in this setting. The dynamics observed here is a subset of the rich behavior expected for soliton-barrier 
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interactions. These dynamics can be extended to more exotic regimes, e.g., by the addition of a species selective 
barrier. The latter will be examined in more detail in our theoretical investigation below. 



III. NUMERICAL INVESTIGATIONS 



In our numerical simulations below, we will assume that the two-component BEC under consideration consists 
of two different hyperfine states of *^Rb, namely the states 1) and |2,— 2) used in the experiment presented 
in the previous section (see also Refs. [TSH21])- In this case, the scattering lengths take the values an — 100. 4ao, 
ai2 = 98.98ao and 022 = 98.98ao (where ag is the Bohr radius); accordingly, the normalized nonlinearity coefhcients 
in Eqs. take equal values: 512 = 322 ~ 0.986. Furthermore, we will assume that the trap frequencies are 

uj± = 27r X 116 Hz and cj^ = 27r x 1.3 Hz, i.e., 51 « 0.0112, and the numbers of the atoms in each component are 
Nd — 70, 000 and Nb = 1, 000 resulting in a chemical potential of approx. 305 Hz for the total atom number. These 
values are similar to the respective ones used in experiments [T9H24] . 

Concerning the parameters of the Gaussian impurity potential, the values of Ei, ,j are taken in the interval 
[—200,200], and we fix the value = £6 = e = 3 fim. Notice that we focus here on relatively narrow impuri- 
ties, as for those we have explored the steady state problem [T7] and, as will be seen below, they already present a 
rich phenomenology. An examination of the effect of the width of the impurity will be deferred to a future study. Our 
principal aim in what follows is to study the scattering of solitons at the impurity potential. To do so, we displace the 
solitons from the trap center, using the initial position value xq = —40 /im, which is sufRciently far from a; = 0, so that 
the solitons do not overlap with the impurity. We then "release" the solitons and observe their subsequent interaction 
with the Gaussian barrier, measuring the fraction and observing also the shape of the atoms that are transmitted, 
reflected, and trapped at the potential. We will study, at first, the case Et, = E^, as per our experimental results (cf. 
subsection A below) and then the case where the impurity acts only on one component, i.e., either Ed = 0, Et ^ 
or Ed ^ 0, Eb = (cf. subsection B below). 

At this point, it is relevant to present a sketch of a state diagram in the parameter space {Eb,Ed), as depicted 
in Fig. [2j The different regimes that appear after the collisions are illustrated by colors and they will be discussed 
below. The capital letters A and B correspond, respectively, to a small and large value of the parameters in each 
region; the cases I {Ed = Ei, = E), 11 {Ed = 0, Et, > 0) and HI {Ed = 0, Ei, < 0) correspond to the principal cases 
that we will examine in what follows. For each A and B we will illustrate the contour plots of the densities of both 
components. Furthermore, in the case where Ed — E(, — E where as we will see below the "particle-like" picture 
is most relevant, we will also display an effective potential energy landscape encountered by the DB solitons. This 
amounts to computing the turning point, say xi, of each of our initializations of the DB soliton at position a;o (whose 
potential energy Vh{xo) in a harmonic oscillator trap we can evaluate). Then, (xi, Vh{xo)) is identified as a point in 
the effective potential energy surface. 



A. Scattering of DB solitons from identical impurities 

We start with the case where both impurities are identical, i.e., Ed = E^ = E, either repulsive {E > 0) or attractive 
{E < 0). 

In the repulsive case oi Ed = E^ = E > 0, the simulations reveal the existence of two different regimes. For 
sufficiently small values of the repulsive barrier E, the DB solitons are transmitted (transmission regime, blue color 
in region I of Fig. [2|. For values of E bigger than a critical value, i.e., i? > 20 Hz (for solitons launched from 
xq — 40 fim) the DB solitons are reflected (reflection regime, red color in the same region I of Fig. [2]). Thus, in this 
case, the solitons behave as classical particles: if they have potential energy (recall that the solitons start with zero 
kinetic energy) larger than the height of the barrier then they are transmitted through it, while they are reflected in 
the opposite case. 

This particle-like behavior can be quantitatively described as follows. First, we consider the maximum of the 
effective barrier height that the DB-soliton particle encounters at the defect, which can be defined as: 

AU = Udef-Ua, (7) 

with Udci (Uq) being the normalized potential energy of the DB soliton at the trap center with(out) the defect: 

^ JZdxV,{x)W{xrdx j:°^dxV,{x)W{xrdx 
C^dx\u'{xW ^ nd^\v'{xW ' 
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FIG. 2: (Color online) A state diagram in the parameter space {Ei,, Ed)- The different regimes that appear after the collisions are 
illustrated by colors: (Blue) Transmission regime. (Red) Reflection regime. (Green) Transmission-reflection regime. (Brown) 
Trapping-transmission-reflection regime. A and B correspond, respectively, to a small and large value of the parameters in 
each region: I {Ea = Eb = E),ll {Ed = 0, Eb > 0) and III {Ed = 0,Eb <0). 



In the above expressions, {u' (x) , v' {x)} and {uo{x) , vo{x)} denote, respectively, the DB solitons at the trap center 
with and without the defect, as found numerically by means of a fixed point algorithm, using as an initial guess 
Eqs. ([5])-([6]). Notice that the solutions with the prime are found by keeping fixed the chemical potentials 111^2 of each 
component to those of the solution without defect. 

Figure [S] shows that the critical value separating the reflection and transmission regime (above and below the solid 
line in the figure) is very accurately obtained from the comparison of the potential energy of the DB-soliton particle 
with AU. Note that in the simulations (see red circles in the figure) we fix the initial soliton location xq and vary 
Ed = Ef, in order to determine the critical value which separates the reflection and transmission regimes for this value 
of Xq or, equivalently, of the potential energy of the soliton, i.e., (1/2)0^X0. 

The contour plots of the densities of the dark and bright components corresponding to relatively small and large 
values of i? > (corresponding to I A and IB in Fig. [2]) are illustrated in Fig. |4] (columns I A and IB). The top row 
shows the initial density profiles of the corresponding dark components. It is important to highlight here (as it will 
become also relevant for other cases) that the dark component sustains an increasing density dip as E increases in 
positive values, while it will correspondingly feature a density bump in the case of increasing negative such values. 
This is a feature of the ground state in the presence of the defect, as the latter repels for i? > and attracts for < 
the atoms in the neighborhood of a; = 0. It is clear also by plotting the effective potentials (in the bottom panel of 
the figure) that the DB faces a weak barrier in the former case and its potential energy is sufficient to overcome it. 
On the other hand, the barrier is considerably higher in the right panel, inducing the reflection of the solitary wave. 

We complete this section by noting that in the case where both impurities are attractive, i.e., Ed = Eh = E < 0, 
the DB solitons are always transmitted after the collision, for every value of E, hence there is only a transmission 
regime depicted by blue color in Fig. [2j 



J°°^ dxVo{x)\uoix)\'^dx ^ dxVo{x)\vo{x)\'^dx 



(9) 



B. Scattering of DB solitons from an impurity in one component 



While the dynamical evolution of the case where the defect acts on both components was found to be fairly 
straightforward, the case where the barrier is imposed selectively on only one of the components was found to be 
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FIG. 3: (Color online) A comparison of the effective potential energy barrier ISU encountered by the DB soliton (solid line) with 
the critical energy of the soliton that separates the transmission from the reflection regime, for different values of Ed ~ Eb = E 
(red circles). The very good agreement (especially for smaller values E) lends support to the consideration of this case as an 
example of a classical (solitonic) particle scattering from a barrier. 

considerably more complex. We considered both the case where E'^ = and E), ^ and that where Ej^ and 
El, — 0. The results show that the following two subcases are equivalent: 

a) {Ed ^0,Eb>0) = {Ed <0,Eb^ 0). 

b) {Ed = 0,Eb<0) = {Ed > 0,Eb = 0). 

The case where Ed — and Ei, < represents the existence of an attractive barrier in the bright component and 
absence of impurity in the dark component. The above equivalence can be most easily qualitatively appreciated in 
that case b) , hence we present it for that setting. In particular, when an impurity attractively affects the atoms in the 
bright component, then it favors the "collection" of atoms near the origin. This, in turn, builds a population of atoms 
in that neighborhood which, through the term proportional to 1712 in the equations of motion, provides a repulsive 
barrier for the dark component. Hence, the existence of an attractive well solely in the bright component becomes 
tantamount to having a repulsive barrier in the dark component. An analogous argument can be used to showcase 
that a repulsive barrier in the bright component, through favoring the absence of atoms in its vicinity, creates an 
effective well for the dark component atoms. The above feature is directly evident in the diagram of Fig. [2] hence we 
only focus on each of the representatives of the cases a) and b) above. 

After the collisions, for Ed — 0,Ei, > 0, part of the energy is transmitted and part of it is reflected. We denote 
this as a transmission-reflection regime (green color in Fig. [2| . For small values of Ef, the DB solitons are mainly 
transmitted, and when Ei, is high enough they are mostly reflected. Equivalent results, when the well depth \Ed\ 
increases, are obtained for the case where Ed < and Eh — 0, which represents the existence of an attractive well in 
the dark component and the absence of impurity in the bright one. For small Ei, > 0, this dynamics can be understood 
as follows. As discussed above, the repulsion of bright atoms produces an effective attraction of dark atoms, hence 
inducing an effective potential well, rather than what was anticipated as a potential barrier. It should be noted here 
that this counter-intuitive effect was quantified in the case of a (5-function potential in Ref. [17] . This effect leads to 
the acceleration of the soliton (with a small back-scatter due to the inelasticity of collision with the defect) visible 
in the left panels of the figure. To complete the discussion of Fig. [5] let us briefly touch upon the right panels of 
the figure. This concerns the case of £'(, < (while Ed = 0, namely case b) above). The corresponding situation 
here, when Eb is small presents a repulsive effect for the dark atoms and as such results in an effective barrier. This 



7 




FIG. 4: (Color online) The case of identical repulsive impurities. Top row: initial density profiles of the dark components for a 
small and large value of -B > 0. Middle rows: contour plots of the densities of the dark components |tt(a;)p and of the densities 
of the bright components |i;(2;)|^. The left panels correspond to a subcritical case, featuring full (nearly full) transmission, while 
the right ones are for a supercritical case with (nearly complete) reflection. Bottom panels: the effective potential encountered 
by the solitary wave. The low height of the barrier in the left column enables transmission, while its increase in the right panel 
induces the reflection. 



prediction is also corroborated by the analytical considerations for the (5-function case of Ref. |17j . This, in turn, leads 
to the reflection dynamics observed in the right panels of Figjsj 

We now turn to the case of large barrier strength in Fig. |6j In this setting, there is a fundamental difference in 
comparison to the case of weak barrier presented previously. This consists of the fact that for small Ef,, the defining 
characteristic in the DB-soliton and defect interaction is the nature of the potential for the dark component (which, 
as we saw, was somewhat counter- intuitively the opposite than the one for the bright component). However, for large 
Eh, the nature of the potential for the bright component becomes important and hence in this case, large positive Eb 
also induces a locally strong repulsive potential for the bright atoms. On the other hand, large negative Ei, creates 
a large attractive potential for the bright atoms. However, the latter tends to favor the trapping of the atoms of the 
bright component together with those of the dark component, leading essentially to the formation of a defect mode, 
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FIG. 5: (Color online) Comparison, for a small value of |_Ei,|, between the case with repulsive bright impurity and the case with 
attractive bright impurity, with Ed = 0. Column IIA: repulsive bright impurity with Ef, = lOHz and Ed — 0. Column IIIA: 
attractive bright impurity with Eb = — lOHz and Ed = 0. Top row: the dark components |M(x,t)p. Bottom row: the bright 
components \v{x,t)\'^. 



alongside a partial reflection of the soliton. These characteristics, namely reflection for Eh large and positive and the 
possibility of trapping, along with reflection for Ef, large and negative, are illustrated in the panels of Fig. [6j 

An example of intermediate values of Ei, and their associated dynamics can be found in Fig. [7j In these examples, 
for positive Eb, the impurity leads to partial transmission and partial reflection, but does not enable the possibility 
of trapping at the defect. The latter possibility is explored for Ei, < 0, whereby there is a fraction of atoms which 
is trapped at the defect, while also a considerable fraction appears to be reflected. We notice that in the case where 
there is no impurity affecting the dark component and where there exists an attractive well in the bright one, as the 
well depth \Eh\ increases, two regimes appear alternatively: a reflection regime and a trapping-transmission-reflection 
regime (red and brown colors, respectively, in region III of Fig. [2]). Indeed, this alternation may be quite complex 
and a characteristic example thereof is presented in Fig. |8] where the fractions of atoms transmitted through the 
defect, reflected from it and trapped in the immediate vicinity of the barrier are measured. Those quantities are 
quantitatively described by the transmission, reflection and trapping coefficients. Notice that the time-dependent 
coefficients are defined as: 

\v{x,t)\^dx j:^\v{x,t)\^dx J^Jv{x,t)\^dx 

" r°° I ^ .M2J ' ■^(^) ~ roo I , ,x|2 , ' -^(0 ~ POO , , ,s|n ■ ■ (10) 

J_oo l^(^'*)r^^ J_oo l^(^'*)r^^ J_oo '^^ 

Figure [s] depicts these coefficients at time f, namely i?', T' and B' , with t' being the time where the bright component 
of the DB soliton reaches its maximum excursion. The figure illustrates that the full dynamics features more compli- 
cated resonant type transmission events, as well as alternating windows of predominantly reflection or predominantly 
trapping. These complex scenarios are beyond the scope of the particle analysis provided herein. We do note, how- 
ever, the apparent similarity of these results with the ones obtained in the case of a single component bright soliton 
which scatters off of a quantum well [39]. In the latter case, the variational analysis was already fairly cumbersome 
even for a (5-function potential, while here it is rendered more elaborate by the presence of two-components and the 
Gaussian form of the barrier. A more detailed analysis, perhaps in the simpler ^-function attractive setting based on 
a two-mode variational ansatz would constitute an interesting problem for further studies. Nevertheless, we should 
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FIG. 6: (Color online) Comparison, for a high value of \Eb\, between the case with repulsive bright impurity and the case with 
attractive bright impurity, with Ed ~ 0. Column IIB: repulsive bright impurity with Ei, — 90Hz and Ed = 0. Column IIIA: 
attractive bright impurity with Eb = — 90Hz and Ed = 0. Top row; the dark components |it(a:,t)|^. Bottom row: the bright 
components \v{x,t)\'^. 
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FIG. 7: (Color online) Similar to the previous two figures, but for an intermediate value of \Eb\. 
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point out that this behavior is similar to that obtained for the case where there exists a repulsive barrier in the dark 
component and there is no impurity affecting the bright component. 
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FIG. 8: (Top) Representation of T' (full blue line, top panel) and R' (dashed red line, top panel) and B' (full green line, bottom 
panel) versus Et for the case Ed = 0. 



IV. CONCLUSIONS AND FUTURE CHALLENGES 

In the present work, we have presented a survey of the collisions of dark-bright (DB) solitons with defects. Initially 
motivated by the potential of experimental studies (a prototypical example of which was shown herein), we considered 
the setting of a dark-bright soliton impinging on a defect potential. In the case of a bright soliton hitting a well [311 SO] 
or a barrier [ITJ this theme has been of intense theoretical and even experimental [HI H?] interest recently. 
However, far less has been done in the realm of dark-bright solitons. 

We have shown that in the case of two equal repulsive barriers acting on both components, the DB solitons 
demonstrate a clear classical particle behavior, which involves transmission for weak potentials and reflection for 
strong ones. Similarly, predominantly transmission type events were observed for equal attractive potentials acting 
in both components. 

On the other hand, we illustrated that more complex scenarios can develop in the case where the impurity acts 
only on one of the two components. We categorized these cases, illustrating the analogies of a repulsive barrier in 
the first component with an attractive one in the second component (and vice- versa). We explained the low barrier 
amplitude cases on the basis of somewhat counter-intuitive, cross-component effective potentials and argued that the 
large amplitude cases may be significantly different due to the role of the defect in both components. We showcased 
the complexity of the latter by means of cases containing transmission and reflection, or trapping, transmission and 
reflection together and by monitoring the dependence of the different fractions (of trapping, transmission or reflection) , 
as a function of the barrier amplitude. 

It would certainly be interesting to extend this chiefly numerical (but also experimental) study further. On the 
experimental side, it would be extremely interesting, although more challenging, to engineer potentials that are 
selective to particular hyperfine states, so that some of the predictions proposed herein could be tested. From a 
theoretical perspective, it would be very relevant to attempt to distill a simple setting (e.g. a ^-function potential) 
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where a theoretical study of the above reported phenomenology could be appreciated in more quantitative terms. 
Numerically, it may also be quite significant to appreciate the effect of the width of the barrier, as here we have 
concentrated on the sign and magnitude (and inter-component interplay of the) barrier. Natural extensions may also 
concern the possibility of scattering in higher-dimensional settings and evaluation of the role of transverse degrees of 
freedom therein. 
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